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The analytic relation between the MS and the pole quark 
masses is computed to O(ctg) in QCD. Using this exact re- 
sult, the accuracy of the large j3o approximation is critically 
examined and the implications of the obtained relation for 
semileptonic B decays are discussed. 



I. INTRODUCTION 

Quark masses, being input parameters of the Stan- 
dard Model Lagrangian, play a very important role in 
high energy physics phenomenology. It is believed that 
the quark masses are generated through a spontaneous 
symmetry breaking mechanism caused by a non-zero vac- 
uum expectation value of a Higgs field. Nevertheless, the 
quark masses remain free parameters of the Lagrangian 
and have to be determined by comparing theoretical pre- 
dictions with experimental data. 

One has to keep in mind, however, that quarks have 
not been observed as free particles and the notion of a 
quark mass relies on a theoretical construction. Differ- 
ent definitions of quark masses exist referring to different 
schemes for the renormalization of the QCD Lagrangian 
of the strong interactions. Two prominent mass defini- 
tions are the MS mass m and the pole mass M. The pole 
mass M is the renormalizcd quark mass in the on-shell 
(OS) renormalization scheme, while the MS mass is the 
renormalizcd quark mass in the modified minimal sub- 
traction scheme which is intimately related to the 
use of dimensional regularization J^]. 

These two renormalization schemes are used in differ- 
ent physical situations. In the MS scheme, the renor- 
malized mass and coupling constant depend explicitly on 
the renormalization scale fj, which is often chosen to be of 
the order of the characteristic scale Q of a process. The 
scale evolution of the MS mass and the coupling constant 
can be done accurately by integrating the renormaliza- 
tion group equations, even for relatively low Q. However, 
the MS mass is by construction sensitive to only the short 
distance (Euclidean) aspects of QCD. For processes in 



which the characteristic scale is large compared to the 
quark masses it is therefore advantageous to adopt the 
MS scheme for the quark mass renormalization^]. 

On the other hand, for processes where the character- 
istic scale is set by the mass of a quark in the initial or 
final state, the situation is different as long-distance as- 
pects of QCD become important. Since in such cases the 
quarks are close to the mass shell, the on-shell scheme 
is a natural renormalization scheme and the pole quark 
mass emerges. Explicit perturbative calculations 
have shown that the pole quark mass is an infrared fi- 
nite, gauge invariant quantity^] and for this reason the 
pole mass of a heavy quark has often been considered as 
a meaningful physical parameter with a corresponding 
numerical value. 

It should be anticipated however that the pole mass 
of a quark is not a physical quantity in a truly non- 
perturbative sense since the confiment of quarks in QCD 
implies that there is no pole in the quark propagator be- 
yond the perturbation theory. For this reason it is natural 
to expect that the pole quark mass might be sensitive to 
long distance effects responsible for keeping quarks and 
gluons confined. This expectation was confirmed recently 
when it became clear that the infrared sensitivity of the 
pole mass reveals itself through contributions in the per- 
turbative series for the pole mass that grow factorially 
at higher orders correspondi ng to a singularity close to 
the origin in the Borel plane ]9|jl0[l . This renormalon sin- 
gularity is important for phenomenology since it implies 
that, indeed, the pole quark mass can not be determined 
with an accuracy better than S M ~ Aqcd P,^0|- 

Taken literary, the renormalon analysis applies to 
asymptotically high orders in perturbation theory. How- 
ever, since the location of the leading renormalon pole in 
the Borcl plane can be correctly obtained in the so-called 
large /3o approximation and since in the lowest non-trivial 
order of perturbation theory the terms proportional to (3q 
dominate, it was conjectured that the leading renormalon 
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1 This is demonstrated for instance in the extraction of the 
MS 6-quark mass renormalized at the scale Mz, from three 
jet events at LEP with a relatively small uncertainty [^J. 

2 This result is proven to all orders in the coupling constant 
in Refs. IF 
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dominates the growth of the perturbative coefficients al- 
ready in low orders for the pole mass. 

It is fair to say that this conjecture is an important in- 
gredient in our understanding of theoretical uncertainties 
for various low-scale physical processes such as semilcp- 
tonic decays of B mesons, determination of the b quark 
mass from T mesons spectrum or the top quark thresh- 
old production cross section. In all these cases eliminat- 
ing the pole quark mass from theoretical formulas is an 
important step in improving the convergence of the per- 
turbative series. Nevertheless, it should be kept in mind 
that, up to date, this physical picture has been checked 
against only one non-trivial order of perturbation theory. 

In spite of its infrared sensitivity, the pole mass re- 
mains to be important for processes where the on-shell 
mass definition has technical advantages. The pole mass 
can afterwards be eliminated in favor of a mass that is 
free from long distance ambiguities such as e.g. the MS 
mass. It often happens that the MS mass is not the 
best choice to parameterize the low-scale processes and 
recently other short-distance low-scale masses have been 
proposed [|il"|-|l3]. An attractive feature of these masses 
is that they can be accurately determined from the anal- 
ysis of the low energy data [p^|-^6|. However, for techni- 
cal reasons, the relation of these masses to the MS mass 
goes via the pole mass and therefore a determination of 
the MS mass from the low-scale short-distance masses re- 
quires the knowledge of the conversion formula between 
MS and the pole mass. 

Hence, the relation between the MS mass and the pole 
quark mass is of prime phenomenological importance. In 
the one-loop order this mass relation was obtained in || . 
The first analytical calculation of the two-loop contribu- 
tion was performed in j(| and this result was confirmed 
in Rl. The three- loop coefficient was first estimated in 



Rcf. 17 using the large 0o approximation. Recently a 
more accurate result was obtained using asymptotic ex- 
pansions and Pade improvements fl8|| . 

In this letter we present the analytical calculation of 
the three-loop relation between the MS and the pole 
mass. Using this exact result, we critically examine the 
validity of the large f3 approximation and discuss some 
implications of the obtained relation for semileptonic B 
decays. A more detailed description of our calculation 
will be published elsewhere. 



II. THE THREE-LOOP RELATION BETWEEN 
THE MS AND THE POLE MASSES 

The quark masses in both, the MS and the on-shell, 
renormalization schemes are renormalized multiplica- 
tively and the connection between renormalized and bare 
quark masses is defined as 



m = Z™ S M, 
m = Z° S M, 



(1) 

(2) 



where mo is the bare (unrenormalized) quark mass, and 
Z^ s and Z® s are the renormalization factors for the 
quark mass in, respectively, the MS and the on-shell 
schemes. The relation between the pole quark mass and 
the MS mass is then expressed as the ratio 
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One sees that in order to calculate the relation between 
the MS and the pole quark mass in the three-loop order 
one needs to calculate the mass renormalization factors 
in both the MS and the on-shell schemes in the three- 
loop order. Fortunately, due to the very simple form 
of renormalization factors in the MS scheme, the mass 
renormalization factor Z^ s is presently known to a suf- 
ficiently high order in perturbation theory and can be 
taken from the literature ITsfl. One obtains: 



Z™ s = l 
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where Nf denotes the number of different fermion fla- 
vors. The dimensional regularization parameter e is de- 
fined through e = (4 — D)/2 with D being the space-time 
dimension. 

In comparison to the MS-scheme, the calculation of 
renormalization factors in the on-shell scheme is much 
more involved and for this reason the renormalization 
factor Z^ is known only in the two-loop approximation. 
In the present work we calculate Z® s in the three-loop 
order of QCD, thereby obtaining the MS pole mass rela- 
tion in the three-loop order. 

Let us explain how the on-shell renormalization factor 
Z® s is computed. This renormalization constant follows 
from considering the one particle irreducible quark self- 
energy parameterized as 

t{p, M) = MY,i(p 2 , M) + {p- M )E 2 (p 2 , M ). (6) 

Here M is the pole-quark-mass and it is understood that 
the mass renormalization is performed in the on-shell 
scheme which means that on-shell mass countcrtcrms 
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proportional to MEi(M 2 , M) are inserted in the dia- 
grams of lower order in a. For the renormalization of 
the strong coupling constant we adopt the minimal sub- 
traction scheme 

^ = -~^(-) 2 +(^~^)(-Y + 0(a% (7) 

7T 7T £ \7T/ \ £ Z 2s J \7T/ 

where /3 = 11/4 - 1/67V/ and fa = 51/8 - 19/24JV) 
are the first two coefficients of the QCD beta function. 
To make the connection with the formal on-shell mass 
renormalization factor Z® s we note that the perturbative 
quark propagator can be written as 



Sf(p) 



p-m + £(p, M) 



(8) 



Since the pole mass of the quark M corresponds to the 
position of the pole of the quark propagator one obtains 



Z 



"*--l+Zi(p 2 ,M)\ 



(9) 



which agrees with the form of the on-shell mass countert- 
erms discussed above. 
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Equation <M) provides the simplest formula for the cal- 
culation of Z~ s . One computes £i(M 2 ,M) to the re- 
quired order in perturbation theory taking into account 
the lower order difference between mo and M by calcu- 
lating lower order diagrams with the appropriate mass 
counterterm insertions. Therefore, in order to obtain the 
three-loop on-shcll renormalization factor Z® s one has to 
compute on-shell propagator-type diagrams up to three 
loops. 

The most efficient way to evaluate those multiloop inte- 
grals is to utilize integration-by-parts identities within di- 
mensional regularization. There are eleven basic topolo- 
gies that should be considered^]. They are shown in Fig.l. 
For each of the topologies one writes down a system of 
recurrence relations based on integration-by-parts iden- 
tities |3|p^|. Solving this system, it is possible to show 
that any integral which belongs to the above topologies 
can be expressed through 18 master integrals. Fortu- 
nately, most of these integrals have been computed in 
the course of the analytical calculation of the electron 
anomalous magnetic moment p3fl and can be taken from 
there. As compared to that reference, we need one addi- 
tional master integral and we also need one of the master 
integrals of Ref. |23| to a higher order in the regular- 
ization parameter e. For completeness, we present here 
the results for these two master integrals. Let us intro- 
duce D x = k\, D 2 = k 2 , D 3 = fcf , At = (fci - k 2 ) 2 ,D 5 = 



(k 3 -k 2 ) 2 ,D 6 = k(+2pk 1 ,D 7 = ki + 2pk 2 ,D t 
with p 2 = — 1. Then: 
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D\D 3 DiD b D 6 D 7 Df. 



= C(e) (2^ 2 C 3 - 5Cs) 



where C(e) = [7r( 2 - £ 'r(l + e)] and ( k = CO) denotes 
the Riemann zeta function. We also give the result for 
Iio (we use notations as in 0)) valid to 0(e 2 ): 
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There are two principal checks on our solution of the 
recurrence relations. First, we have computed the three- 
loop anomalous magnetic moment of the electron and 
confirmed the result of Ref. |^|. Second, the actual cal- 
culation of the on-shell quark mass renormalization con- 
stant Z® has been performed in an arbitrary covariant 
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FIG. 1. Examples of three-loop quark propagator diagrams 
corresponding to eleven integration topologies. 



For Abelian gauge theory, like QED, the number of basic 
topologies would be just four. 
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gauge for the gluon field. The explicit cancellation of the 
gauge parameter in our result for Z^ s is an important 
check of the correctness of the calculation. 

Having computed the three-loop contribution to Z^ s 
in this way and using Eqs.([||]) we obtain the three-loop 
relation between the MS quark mass and the pole mass. 
Below we present the exact result only for this mass rela- 
tion but we should note for future reference that Z® s is 
easily recovered using Eq.(j3[). The relation between the 
masses is expressed in terms of color and flavor factors 
as 



m{M)= M 



1-Cf 



+c A 4 2) 



t r n l 4 2) 



(t)'(< 

T r N h 4 2) 



4 2) 



+C F T R N L df - 
+C A T R N H df ] 



(3) 
F "1 



Cld 



CpC A <^2 



(3) 



(3) 



C F T R N H df - 
- T 2 R N L N H 4 3) 



+T R Nld$ 



o[2- 

7T 



T C A u 3 

C A T R N L df 
(10) 



where Cf and C A are the Casimir operators of the funda- 
mental and the adjoint representation of the color gauge 
group (the group is SU(3) for QCD) and T R is the trace 
normalization of the fundamental representation. Nl is 
the number of massless quark flavors and Nh is the num- 
ber of quark flavors with a pole mass equal to M. Also 
a s = a\ NL+NH \M) is the MS strong coupling constant 
renormalized at the scale of the pole mass /i = M in the 
theory with Nl + Nh active flavors 4 . 

Our calculation yields the following result for the co- 
efficients 4 n) in Eq. (pi): 
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where a 4 = £ l/(2"n 4 ) = Li 4 (l/2) w 0.517479. 

n=l 

For standard values of the QCD color factors Cf = 
4/3, = 3, T R = 1/2, and assuming the number of 
heavy flavors Ajy to be equal to one, the result reads: 
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4 We should note that, if needed, Eq.([lo|) is easily expressed 
in terms of m(/i), a s {fi) and log(/i/M) by using the renormal- 
ization group equations for m(fj,) and a s (/j). 



Numerically, we obtain: 



m(M) 



4 /a 



M =1 -3 1VJ + IT) d-0414 AT, -14.3323) 



4 



+ (— ) 3 (-0.65269 Nl + 26.9239 N L - 198.7068) . (12) 

We note that the exact three-loop 0{N L ) and 0(iV°) 
coefficients are within the errors of the values 27.3(7)7Vl 
and 202(5) obtained in ||l8fl, nevertheless for the indi- 
vidual color structures CpC\ and C\Nl the difference 
between the exact result and the numerical values in |Q 
is larger than the error bars quoted in that reference. 

From Eq. flll|) a relation between the pole mass M and 
the m(m) , the MS mass normalized at the scale fi = rn(fi) 
is easily derived. One obtains 



m(m) 



4 /ots 

3 V 7T 



(—J (-!-0414 N L + 13.4434) 
(— ) J (0.6527 Nl - 26.655 N L + 190.595) , (13) 



with a s = a s (rn). 

It is interesting to look at these results when they are 
parameterized in a way, that separates the effects related 
to the running of the coupling constant (i.e. contribu- 
tions that are related to a change of the scale /j, in the 
coupling constant) and the remaining part. To do this in 
a consistent manner, we first use the decoupling relations 
to reexpress the MS coupling constant a s defined in the 
theory with Nl massless +1 massive flavors through the 
coupling, defined in the theory with Nl massless quarks 
only pg , Bl| . Using the first two coefficients of the /?- 
function: 



f3 1 = ±( W2 -™NL 
h 16 V 3 



we obtain: 
M 



1 + ^(v) + (v)' (6 - 248 /3o_ 3 - 739) 



m(m) 



+ (— ) (23.497 Pi + 6.248 (3 X + 1.019 O - 29.94), (14) 

where a s now stands for the MS coupling defined in the 
theory with Nl fermions only. 

From the above equation one sees that the large /3q ap- 
proximation works well at O(af) especially for Nl = 3, 
but in part because of a cancellation between the confor- 
mal (i.e. /3-independent) term and the term proportional 
to 0i. The magnitude of conformal terms is large and 
they exhibit a very fast growth from one order of pertur- 
bation theory to the other. 

In order to check that the magnitude of conformal 
terms is not an artifact of our use of the MS coupling 
constant as an expansion parameter, one can rewrite the 
equation for the mass ratio in terms of the ^-scheme 
coupling ay = av(m). To the necessary order the rela- 
tion between ay and a s is given in Ref. |p4| . Rewriting 
Eq.(pT[) in ay, one obtains the result which does not 
improve as compared to the MS scheme. The overall 
magnitude of corrections gets somewhat decreased, but 
conformal coefficients remain very large and grow rapidly. 



This situation makes it difficult to interpret the high 
accuracy of the large /3q approximation at 0(ot s ) simply 
from the assumption of a leading infrared renormalon 
dominance^]. On the one hand we clearly see that the 
large f3$ approximation remains extremely successful in 
the third order of the perturbative expansion for the pole 
mass. On the other hand, part of the reason for this 
success seems to be a cancellation between subleading 
/3-dependent and fast growing conformal terms. Since a 
different physics is usually associated with these contri- 
butions it is unclear how this situation will extrapolate to 
the next order in the MS pole mass relation. We refrain 
from drawing a definite conclusion on this point at the 
moment. 

Since exact calculations in even higher orders of per- 
turbation theory seem formidably complicated, it is per- 
haps reasonable to look for other approximation schemes 
which, as compared to the large /3q approximation, will 
preserve in a better way the field theoretical structure 
of QCD. One possibility is the large N c approximation 
|p5| . Its attractive feature is that only planar diagrams 
contribute to leading order in N c and this significantly re- 
duces the technical burden of the calculations. Since in 
realistic QCD we are interested in the applications where 
N c ~ Nl, one can go a step further and consider a situa- 
tion where both the number of colors N c and the number 
of light fermions Nl are considered to be large numbers. 
The correction to this approximation is expected to be of 
the order 0(N~ 2 ) which for iV c = 3 is about 10 percent^. 
The result for the MS mass in this approximation reads: 



m(M) 
M 

-m 



= 1-1.5 



£) + £)W*-ie.is) 

(15) 



-0.73 N[ + 30.41 N L - 222 



Comparing this equation with Eq. (|lj) one sees that the 
large N c approximation works with the expected accu- 
racy and hence might serve as a useful starting point for 
the analysis at an even higher order. 



We note that the large j3o limit and the infrared renormalon 
dominance are not the same thing. Also in the renormalon 
context corrections to the large /3q limit grow with the order 
in perturbation theory and eventually must make the large /3o 
limit unreliable. We thank M. Beneke for correspondence on 
this point. 

6 For contributions of massive fermion loops the accuracy of 
the approximation is formally 0(N~ 1 ) but due to the mass 
of these fermions these contributions are usually small. 
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III. THE THREE-LOOP ANALYSIS OF 
SEMILEPTONIC B DECAYS 



Let us now turn to the analysis of the semileptonic de- 
cay process B — > X u eu e . We want to investigate whether 
various existing estimates of the 0(a\) correction to this 
process are compatible with the notion of an improved 
quality of the perturbative series at low orders when the 
pole quark mass is eliminated in favor of a properly de- 
fined low-scale short-distance mass. Some of these esti- 
mates refer explicitly to the MS mass and this investiga- 
tion requires the use of the three-loop quark mass relation 
Eq. @. 

Recall, that the Heavy Quark Expansion demonstrates 
that the non-perturbative effects are small in this pro- 
cess [^6| and an important issue is the accurate calcu- 
lation of the perturbative corrections to the quark level 
decay b — > uev e . To 0{a 2 a ) such a calculation has been 
performed in Ref. |27[| . Recently the three- loop correc- 
tion to that process has been estimated in Ref. |28| using 
asymptotic Pade approximants. It is unclear to us at the 
moment how accurate that estimate is, but we will take 
the number quoted in Ref. |^8| as a reference point. For 
five active flavors we write: 



r sl = 



G 2 F \V ub \ 2 m(m) 




-(5)' 



r,(5) ' 



1 + 4.25 



26.78 



(16) 



The central value of the third order coefficient is due to 
Ref. |p8| and A parameterizes possible uncertainties in 
this estimate. For the chosen normalization point fi = to, 
one sees a rapid increase in the coefficients of the series 
in Eq.(^6|). The reason for that is thought to be a poor 
choice of the renormalization scale, since a typical energy 
release in B decays is of the order m b /5 rather than m b . 
Since at very low scales the logarithmic running of the 
quark mass is considered unphysical, it was suggested 
to use other short-distance masses which have a more 
transparent meaning at a low normalization point. In 
this paper we would like to check the behavior of pertur- 
bation series once the decay width is expressed in terms 
of the so-called 15 mass. Strictly speaking, the 15 mass 
introduced in Ref. Q is not a truly short-distance quan- 
tity and its usefulness for B decays when both pertur- 
bative and non-perturbative corrections are taken into 
account remains to be proven. Nevertheless, we use it 
here because its relation to the pole mass is known to a 
sufficiently high order in the perturbative expansion. It 
was suggested that the correct way to incorporate the 15 
mass into semileptonic B decays is the T expansion |29[| 
and we will denote different terms in this expansion by 
a formal parameter e (which is eventually put to one). 
We use the relation between the MS mass and the pole 
mass to rewrite Eq.([l^) in terms of the 15 mass (via its 



relation to the pole mass) and check if the result is con- 
sistent with an improved behavior of the perturbation 
series. Below we consider four possible scenarios. 

i) The large Bq estimate of the three-loop coefficient of 
r s i is accurate^. Using the mass relation, derived in this 
paper, we find that this requires A w —100 in Eq.(fl6l). 
We then rewrite the decay rate in terms of the 15 mass 



and obtain for a 



(4) 



r H , = 



G 2 F \V ub \ 2 m ls 5 

1927T 3 



0.22: 



1 - 0.115e - O.QSle 2 - 0.027e 3 



ii) Let us assume that the pattern of the n-th loop 
coefficient in Eq. (jll^) is well approximate d by 5™ (see Ref. 
filf ) . For the three-loop coefficient in Eq. ([lq ) this implies 
A « —75. Expressed in terms of the pole mass the three- 
loop coefficient of the width is then a factor 0.9 smaller 
than the large Bo estimate. Expressed through the 15 



mass, the decay width for a 



(4) 



0.22 reads: 



r = G 2 F \V ub \ 2 m ls 5 . _ Q _ Q 031e2 _ Q 0i8e3 _ 

1927T J v ' 

iii) Assume that the estimate of Ref. j28| is accurate. 
Expressed in terms of the pole mass the three-loop co- 
efficient of the width is then a factor 0.7 smaller than 
the large Bo approximation. Rewriting the decay width 
through the 15 mass and using cSf^ — 0.22 we obtain: 



r a , = 



G 2 F \V ub \ 2 m ls * 

1927T 3 



(1 - 0.115e - 0.031e 2 + 0.008e 3 



iv) One can see that at the two-loop level the exact 
result for the decay width (in terms of the pole mass) 
is about a factor 0.8 smaller than the large Bq approxi- 
mation. Let us imagine that in the three-loop order the 
exact result is about half the size of the large Bo approx- 
imation. We note that this translates into A « 80 in 



Eq.fllJ). For aT> = 0.22 we obtain: 



r a , = 



G 2 F \V ub \ 2 m ls 5 
192vr 3 



1 - 0.115e - O.OSle 2 + 0.035e J 



We see therefore that there is a window for the pa- 
rameter A where a parameterization of the decay width 
through the 15 mass seems to succeed in making the 
perturbative series convergent. However, there is one un- 
pleasant feature. The magnitude of the 0(e 3 ) term in the 
T expansion strongly depends on the value of the strong 
coupling constant and this is because various terms which 
scale differently in a s contribute to a given term in the T 



7 When we refer to the results of "the large f3o approxima- 
tion" in this context, we have in mind the expression for r a i 
written in terms of the pole mass. The predictions for F s i in 
the large Bo approximation can be found in Ref. [|30(] . 
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expansion^. In general, this dependence is quite strong. 
For example, by using a s = 0.25 we would have obtained 
0.06 as the coefficient of the e 3 term in the fourth exam- 
ple. 



IV. CONCLUSIONS 

The main result of this letter, an analytic three-loop 
relation between the MS and the pole masses, is given 
by Eq. (|TT]) . It provides a first step towards a more pre- 
cise analysis of QCD effects in different low-scale shortx- 
distance processes. 

Although this result is in good agreement with the 
large @q approximation, we have argued that it is difficult 
to understand the structure of subleading terms. Never- 
theless, our analysis of semileptonic B decay B — > X u ev 
shows that it is likely that the general physical picture 
of improving a badly behaved perturbative series at low 
orders by eliminating the pole quark mass remains valid 
also in the three- loop order. However, to answer this 
question unambiguously, one would have to perform a 
three-loop QCD calculation of semileptonic decays of a 
heavy quark and this seems to be a very non-trivial task 
at present. It may be that the large N c approximation 
will provide enough technical simplifications to be feasi- 
ble beyond the presently known perturbative orders and 
at the same time be accurate enough to be meaningful. 
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